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Abstract 

Numerical studies of the Anderson transition are based on finite-size scaling analysis 
of the smallest positive Lyapunov exponent. We prove numerically that the same scaling 
holds also for higher Lyapunov exponents. This scaling supports the hypothesis of the 
one-parameter scaling of the conductance distribution. From collected numerical data for 
quasi one dimensional systems up to system size 24? x oo we found the critical disorder 
16.50 < Wc < 16.53 and the critical exponent 1.50 <v< 1.54. Finite-size effects and the 
role of irrelevant scaling parameters are discussed. 

PACS numbers: 71.30.+h, 71.23.-k, 72.15.Rn 
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The main problem of the theory of Anderson transition is to prove that there is only one 
relevant parameter which controls the behavior of all quantities of interest in the neighborhood 
of the critical point. An excellent example of such a quantity is the smallest positive Lyapunov 
exponent (LE) zi which follows the one-parametric scaling relation 

z,{L,W) = f{L/aW))\ (1) 

In (H]), is the disorder, L defines the width of the quasi-one dimensional system LxLxL^ and 
^iyV) is the universal scaling parameter. In QID geometry, Lyapunov exponents Zi are defined 
through eigenvalues ti of the transfer matrix T^T as Zi = 2-^logti. In the limit >> L, all 
zs are self-averaged quantities [0. Relation (|^) determines the disorder and the system size 
dependence of zi in the neighborhood of the critical point Wc and enables us to determine Wc 
and critical exponents for conductance (s, W < Wc) and for localization length (z/, W > Wc) 

[0 i i i i 

In the pioneering work critical parameters were found as Wc ~ 16.5 ± 0.5 for the box 
distribution of the disorder energies, and s = u ^ 1.5 ± 0.1. This result was later confirmed 
by more accurate numerical studies 0, H, ^, and also by analysis of the level statistics 0. 
Calculations performed for different microscopic models confirmed the universality of exponent 
u within a given universality class . 

To complete the proof of the universality of the metal-insulator transition, the one parameter 



scaling should be found also for more realistic variables, such as the conductance g [|1^, [TJ]. This 

must be done for the cubic samples. Here, owing to the absence of self-averaging, it is necessary 

to test the universal scaling of the whole distribution P{g). It is unrealistic to perform such 
^Instead of zi, the inverse quantity A = 2/zi is commonly used. The present discussion is identical for both 

quantities, but z's seem to be more natural variables for our purposes. 



analysis with the numerical accuracy comparable to that achieved from QID studies. Therefore, 
previous studies of P{g) concentrated only on estimation of the conductance distribution at 
the critical point [|1^, |1^, [1^ and in the metallic and localized regime . 



The aim of this work is to support the idea of the one-parametric scaling of the conductance 
and of its distribution. Instead of the study g, we prove numerically that the higher Lyapunov 
exponents Z2,Z3,... follow the same scaling behavior as the first one in the QID systems. 
Common scaling proves that the spectrum of the transfer matrix in the QID limit is determined 
only by one parameter. Strong correlation of zs gives also the serious basis for the generalization 



of the random - matrix theory to the description of the critical region ||T6|, |17|, . 

Although we deal only with QID geometry, it is reasonable to suppose that the observed 
correlation survives also for the cubic geometry. Then the relation between g and zs, g = 
J2i cosh.~'^ Zi/2, [|19| assures that g also follows the one-parametric scaling. 

Collected numerical data also provide us with a very accurate estimation of the critical 
disorder Wc and the critical exponent i^. It is the first time that numerical data for system size 
L > 16 has been collected and analyzed. Our data for large L enable us to check the finite-size 
corrections to scaling proposed in p. 

The scaling behavior of higher LEs was originally studied in the Henneke's PhD Thesis 
[0 . Due to the insufficient accuracy of his data and small system size, no acceptable proof of 



the common scaling was found. The first indication of the common scaling was shown in 



and generalized to the neighborhood of the critical point in |]16|, ^Tj. The common behavior of 
higher LEs, ~ i is well known in the metallic regime; it was already used in to explain 
the physical meaning of the scaling parameter C,{W), and confirmed by random- matrix studies 



||19|| . In the localized regime, zs follow the relation Zi{W, L) = zi{W, L) + Aj with W— and L- 



independent constants Aj ||T6 . 



For the QID systems x we calculated all LEs for cca 21 different values of disorder, 
16 < ly < 17. L grows from L = 4 up to 24. For the smallest L, the relative accuracy of the 
first LE zi{W, L),ei = y/'^^Efl^/ zi{W, L) was 0.05% while Si was only 1% for L = 21, 22 and 24, 
being 0.5% for L = 16, 18. The accuracy of higher LE is much better; in particular, £2 ~ ^i/2 
and Eg ~ O.lT^i for all system size. 

The interval of the disorder is narrow enough to approximate the W dependence of z^s by 
the linear fit 

z,iW,L)=zf\L) + Wzl'\L), J = 1,2,... (2) 

Small differences between fits containing higher powers of W and (|]) appear only for L > 18 
and even then they do not exceed the numerical inaccuracy of the raw data. The typical W - 
dependence of our data is presented in Fig. 1 for 22 • 



The scaling behavior requires that |20 



Zj{W, L) ^ + Ax {W -W,) X L"", a = 1/v. (3) 

Comparison of (0) with (|^) offers the simplest way to estimate the critical exponent a. In Fig. 
2. we present the L-dependence of z^p for the first six LEs and for zq. It confirms that close 
to the critical point these LEs scale with the same exponent a: 

a « 0.655 ±0.010 (4) 

which determines v = 1.526 ± 0.023. This estimation is in very good agreement with the result 
of MacKinnon and differs slightly from |^ . 
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Figs. 1 and 2 also show also the important influence of the finite-size effects (FSE) in the 
present analysis. Evidently, the small L data are of no use in the analysis of higher LEs. We 
found that the numerical data for Zj could be used only when 

L> Zj. (5) 

It is easy to understand. If Zj > L then the jth channel is rather "localized" than critical 
on this length scale. Therefore only a small part of the spectra which fulfills the relation (^) 
follows the scaling behavior. The rest of the spectrum depends on L even at the critical point. 
This conclusion is supported also by the analysis of the density p{z) of all LEs for the cubic 



samples At the critical point, the number of system-size independent LEs grows as ~ L 



when L — > cxD |jT6|. As zi ~ 3.4, the above mentioned effect does not influence the analysis of 
the first LE zi. Nevertheless, other FSE must be taken into account in the scaling analysis of 



More reliable estimation of the exponent a and of the critical disorder Wc is given by 
the position of the minimum of the function 

F{Wc, «,•••) = ^ E Zy^.l^m L) - zf{W, L)f. (6) 

W,L C^ill^'^J 

In (H), = J2w,L is the number of points, and . . . stays for all other fitting parameters. 

The natural choice of the fitting function zj^ in @ is the rhs of (0). None FSE are explicitly 
included in (^. Nevertheless, it still enables us to test the sensitivity of the critical parameters 
to the size of the analyzed systems. To do so, we considered different sets of input data Zj{L, W) 
with the restriction Lmin < L < L^ax {L^in < 12). Then, the Lmin- and Lmax- dependence of Wc 
and a was analyzed. While the influence of the choice of L^a.x is, as supposed, negligible, both 



Wc and Zjc increase with Lmin- We found the Lmin - independent resuhs only for the two smallest 
LEs Zi and Z2. For the higher LEs, critical parameters do not reach their limiting values even 
for Lmin = 12. In difference to Wc, the estimation of the critical exponent a does not depend 
on the choice of interval of L. Obtained data are in good agreement with estimation (|^) for all 
LEs under consideration. 

The weak Lmin - sensitivity of the critical exponent agrees with an assumption that FSE 
influence primarily the py-independent part of Zj [Q. Fig 1. offers a simple interpretation of 
this result: to eliminate FSE one has to shift each line by the disorder- independent constant 
B{L) which should be added to the rhs. of @. The proper choice of B{L), assures that all 
lines cross at the same point as it is proposed by the scaling theory. Finite size corrections to 
the line slope are only of the " higher order" . 

Slevin and Ohtsuki p fitted zi{W, L) (resp. its inverse zf^) to the more general function 

zf\L, W) = Z,c + E E Anmx'^y'^ (7) 

n=0 m=0 

with = 3, Ny = 1. In (0), x = {w + hiw^ + 62w3)L", w = W - and y = L^^ with 
/5 < 0. Exponent (3 represents the second critical (irrelevant) scaling exponent. We applied 
this function to our data with restriction (|^) and with hi = b2 = 0,n + m< 1. Then 

zf{L, W) = Zjc + Ax{W - Wc)L'^ + BL^. (8) 

We have checked that more sophisticated fits do not provide us with any reasonable improve- 
ment of the accuracy of critical parameters. 

To test the quality of the fit (H), we again studied the sensitivity of our results to a change 
of the input data. Evidently, for large enough Lmin the role of the irrelevant scaling exponent is 
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negligible. The finite size effects become small and difficult to measure. Value of the irrelevant 
critical exponent (3 obtained from fitting function (|) decreases to ~ —20 for large Lmin- 

For small values of L^[^, however, the three-parametric fit (Q) still does not provide us with 
the Lmin-independent estimation of critical parameters. We averaged therefore the values of 
Wc and a obtained from various choices of I/min- 

Table 1. presents our results for the first five LEs obtained from fits @ and (§). On the 
basis of the presented data we estimate 

16.50 <Wc< 16.53 and 1.50 <u< 1.54 (9) 

These values are in a very good agreement with |Q. 

Our results (^) differ from that obtained by many parametric fitting procedure in Ref. 0] 
(Table 1.). None of the analyzed statistical ensemble provides us with such high value of u. 
This discrepancy is probably caused by different input data. Contrary to previous treatments 
1^, we collected data for large system size in order to simplify the fitting function. The main 
shortcoming of this strategy is a lower accuracy of our data for Zi. On the other hand, the fact 
that the results obtained from the many parametric fitting procedure depend on L^in indicates 
that the fitting function (§) is still insufficient to reflect completely the corrections to scaling. 
The only way to obtain a more accurate estimation of the critical parameters is to collect more 
exact numerical data for large system size. 

To conclude, we have collected numerical data for the quasi one dimensional Anderson model 
up to system size L = 24. Our data prove that higher Lyapunov exponents of the transfer matrix 
follow the one-parametric scaling law. The critical exponent u coincides with that calculated 
from the scaling treatment of the smallest LE. The scaling holds only for Lyapunov exponents 

7 



which are smaller than the system size considered. 

The common scaling enables us to express all relevant LEs as the unique function of the 
first one. Evidently, the same holds also for any function of zs. This indicates the validity of 
the scaling theory for the conductance. However, our analysis was restricted to the quasi-one 
dimensional systems. Rigorous proof of the one-parametric scaling of the conductance requires 
repeating the performed scaling analysis for the cubic samples, where no self-averaging of zs 
and of g takes place. 

We show for the first time, that the numerical data for higher LE could be used for cal- 
culation of critical parameters of the metal- insulator transition. The numerical accuracy of 
higher LE is much better than that of zi. The price we pay for it is a stronger influence of the 
finite-size effects which causes that the data obtained for small system size cannot be used for 
the scaling analysis. The best compromise between the accuracy and FSE offer data for the 
second LE Z2- We discussed the methods of elimination of the finite size effects and estimated 
the critical disorder and the critical exponent by relation (j^). 
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j 


^min 








a 


h' 


13 


1 


4-5 


24 


16.515 


3.46 


0.644 


1.55 


-3.5 


1* 


8-12 


20-24 


16.505 (10) 


3.451 (07) 


0.681 (15) 


1.470 (30) 




2 


5-10 


24 


16.527 (02) 


5.588 (02) 


0.654 (08) 


1.529 (18) 


-3.2 (6) 


2* 


10-12 


22-24 


16.500 (07) 


5.500 (07) 


0.659 (05) 


1.517 (11) 




3 


9 


24 


16.508 


7.167 


0.647 


1.545 


-6.0 


4 


8-10 


24 


16.504 (02) 


8.422 (05) 


0.663 (04) 


1.509 (9) 


-3.7 (2) 


5 


9-12 


24 


16.517 (16) 


9.560 (30) 


0.661 (06) 


1.513 (14) 


-3.3 (8) 


li 


4 


12 


16.500 (50) 






1.515 (33) 




li 


6 


12 


16.448 (14) 






1.59 (3) 




li 


4 


14 


16.540 (10) 






1.57 (2) 


-2.8 (5) 


l*i 


8 


14 


16.514 (07) 






1.58 (5) 





Table 1: Critical disorder Wc and critical exponent u as have been found from numerical data 
for the jth LE for the three-parametric and two-parametric fits (*) and their comparison with 
other results. Number of analyzed points is ~ 21 x (Lmax — -^min)- The minimum of F was 
found < 1.05 for all analyzed sets (with exception of z^, where it was 1.09). 
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Figure 2: The L-dependence of zf^ for the first six LE and for (counted from below). The 
slope determines critical exponents as a = Xjv. Inset: values oiXjv found from presented fits. 
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